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In this paper, frequency oscillation has been present as it is an important 
issue in power systems, especially when it is joined to the new trend to use 
alternative energies as an energy source, as well as to be a big risk for the 
inter-connection of the modern electric power networks. Wind farm acts as 
alternative energy and connected to two buses on the Iraqi power system. 
Because the low-frequency oscillation monitoring needs be accurate and 
fast, the main objective is to propose a novel online monitoring system 
consisting of phasor measurement unit's (PMU) with artificial intelligence 
neural network (PMU-NN). The location of the phasor measurement units 
has been optimized using (graph-theoretic procedure algorithm) and the 
function for the artificial intelligence (NN) is radial basis function 
(RBFNN). The data information from phasor measurement units is the 
inputs to the artificial intelligence system then predictions are made 
Information on low-frequency oscillation (target). The MATLAB toolboxes 


(PSAT & NN) used to obtain results. Finally, from the results, the validity of 
the proposed (PMU-NN) system has been proven and tested on the Iraqi 
power grid (24 bus) in several cases and several places on the network and 
the comparison was made with the analysis model. 
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1. INTRODUCTION 

In modern electrical power systems, and with the development taking place in the world of 
industries that require ensuring the equipped capacity of these industries in terms of availability, quality, 
reliability, and protection, it has become imperative that this equipped capacity be monitored by an intelligent 
and rapid response system capable of predicting and giving very accurate data to ensure appropriate decisions 
are made in a failure occurs. The PMU's is most important devices used in electrical power systems network 
through which it is possible to evaluate any site for any event that occurred or will receive the network [1]— 
[3]. Previously, the case was evaluated without the need for a unit of measurement. But now, with the unit of 
measurement, the bus calculated as neighbors, as it makes all observations and measurements for both 
voltages and currents are available [4], [5]. In fact, in the major electrical power systems, all things are 
financially costly, with the installation of phasor measurement unit within the network devices, and because 
of its financial costs in terms of infrastructure and installation costs, and this is why major companies avoid 
installing it in all places [6], [7]. When alternative energies are used within the power system and represented 
by modern energy systems, the reaction of the power networks during event of a small breakdown fluctuating 
as if voltage regulators were used, which will reduce the risk of instability in power system [8], [9]. The main 
step to monitor oscillation is the power system electromechanical identification modes. We can divide the 
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identification of technologies into two groups of methods, the first is modal analysis methods and the second 
is modal identification [10], [11]. Modal analysis methods are one of the greatest approaches for offline 
study, but for time-overwhelming and not performing well for very large systems. Modal identification 
methods of identification comprise analysis, pattern analysis meters, non-linear and non-fixed methods of 
analysis [12]. There is no size restriction nevertheless this requires some time field simulations to excite and 
differentiate all between the region's modes in the network [13]. Now that technology has advanced using 
everything necessary to make the power system a perfect system, therefore, the development of on line 
monitoring system has developed very important and useful to provide this goal [14], [15]. 

Finally, in this paper a new online monitoring system is built, (PMU-NN) proposed to achieve the 
objective of this paper, MATLAB toolboxes (PSAT and neural network) are used and tested on IRAQI 
power system (24 bus). The contributions of this paper as: i) Proposing a linear model of PMU as a valuable 
technique to get the input variables and then feed the neural network system to estimate the information of 
oscillations; ii) It is the first time to use computational methods such as integrated (PMU-NN) on the Iraqi 
power system as an online monitoring system; and iii) Providing a realistic integrated system PMU-NN 
within proficiency to be applied in big scale cases. 


2. FORMULATION OF POWER SYSTEM INDICATORS 

There are many indicators to expect the power systems state such as frequency oscillation, participating 
generation's units, mode index and localness. Power system indicators are measures used to assess the health and 
stability of a power system. Frequency is a critical power system indicator and deviations from the standard 
frequency can indicate potential problems. Participating generation units, mode index, and localness are 
additional indicators that can provide insight into the stability of the power system. These indicators help 
operators and engineers make informed decisions about managing and maintaining the power grid. 


2.1. Damping index 

As is expected in electrical power networks, the fluctuation that may be happens as a consequence 
of any sudden occasion is variable or that it is in varying degrees from one case to another according to the 
strength of the sudden event. The electric power systems monitoring system must be able to determine the 
strength or weakness of the impact value for that sudden event, for this the damping index (DJ) is used [16] 
as in (1). The damping index rang between (0-1) if it was equal to 0 which means the system is very good 
damped if it is more than 0 to => 1 which means well-damped, else it means poorly. 


; “) 
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Where: ¢ :is the ratio of damping; ¢ min,admissipie : 18 the min. acceptable damping limit; n is index norm; and 
Cabs 1S an absolute of the damping limit. 


2.2. Mode index (MI) 

The frequency ranges and the participating generation units, it is the points to show the differences 
between the local/inter-area modes. However, in the power grid system that have a large of buses the 
observation may difficult to define the mode if it is global or local. Local modes from 0.8 to 2 Hz include 
swinging the components at the generator concerning rest of power grid. Inter-area modes from 0.1 to 0.8 Hz 
are correlated with swing of machines in area of the power system against machines in other areas. MI ranks all 
the modes rendering to their localness [17], [18]. If the MI value is high that means it is a local mode or if it has 
a small value that is mean an inter-area mode. The MI of ith electromechanical mode is given in (2). 


MI, = "(1 — py)” (2) 


Where: I,, = Electromechanical mode; Nz = generators No.; pj; = normalized participation factor value; 
and n = clustering property. 


3. PROPOSED METHODOLOGY (METHODS) 
3.1. Optimal placements of phasor measurement units 

Graph-theoretic procedure algorithm (GTP) [19], [20] has been used to reduce PMU devices that 
can connect to the Iraqi power network (24 bus) for full network observable. The details of (the 24 bus) Iraqi 
power system has been shown in Table |. Figure 1 shows the algorithm of optimal placement. Figure 2 
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shows the single line diagram for an Iraqi power system with PMU installation and two wind farms are 
installations on 1-MUSP and 15-BGE4. 


Table 1. Iraqi power system details (24 bus) with PMU bus location 


Name Bus NO LINEs Gen NO Load NO 
Iraqi power 24 36 11 19 
system PMU’s NO. PMU’s Location Time in (Sec.) Requirement 
8 4-BAJG, 10-HRTP, 12-MSL4, 13-BGW4, 0.36341 


16-BGN4, 18-BGC4, 20-KUT4, 21-QIM4 


Place PMU at bus in the 
unobservable area 


Determine Systems 
Converge with PMUs 


Is System 
Observable 
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19-DYI4 
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Load } 
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Figure 2. Iraqi power system single line diagram with PMU installation and wind farms [23] 


3.2. Artificial intelligent neural network 
The neural network acts as a process box where it manages inputs mathematically [24]. As it selects 


the mathematical operations in a way that suits the entry (variables) and works to create parallel relationships 
with all elements. The approximate values produced from these parallel relationships vary with precision. 
These results are obtained by training input variables. The artificial intelligence radial basis function 


A novel online monitoring system of frequency oscillations-based intelligence ... (Husham Idan Hussein) 


1592 O ISSN: 2088-8694 


(RBENN) has an input, hidden, and output layer (target). In this paper, the wide function used is the radial 
basis function (RBF) that produces the target ym (predicting or estimating the LFO information). 


n= oxy (ul 
y p 202) 


(3) 


Where: (X; — U;): square distance between the input feature vector (X;) and center vector of (U p hidden node 
as measured by some norm for that radial basis function. Here in this paper, the hidden layer number is three 
of size 18 in RBFNN. The input vectors define as phase voltage values and angles as shown in Figure 3. 


Hidden 
(RBF) layer 


Input layer 


Figure 3. Radial basis function neural network (RBFNN) [4] 


4. PMU_NN INTELLIGENT MONITORING SYSTEM 

In this part of the paper, the phasor measurement units provide the input data, which includes both 
voltages and their angles [25], while the artificial intelligence (NN) works to treat the inputs as variables 
multiplied by weights and chooses a mathematical method and layers to extract the targets or the outputs and 
it is more than one output such as (MI, DR, and FI) mode index, damping ratio and frequency respectively, 


for each iy, mode, as shown in the Figure 4. Figure 5 shows the general algorithm for the proposed work 
(PMU-NN) based monitoring system. 


Artificial Inteligent Output 
Machine (Target) 
= 


¥ 
7 ‘N 


Data Generation from 
PMU in to NN 


Power System Operations 
Randomly Loading (i) 


Divide the Data set for 
(MI, DR and fi) 


All 
Operation 
Condition 
Done 


Training and Testing 


Data for Purpose 


Figure 4. Block diagram for PMU-NN Figure 5. Flow chart of PMU-NN algorithm-based 
monitoring system monitoring system 
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5. RESULTS AND DISCUSSION 

The innovative technique in this paper has been tested on Iraqi network with (24 buses), 
(11 generators) and (17 lines) with the following features: The Iraqi power system high voltage 400 kV, 
popular base power of 100 MVA. A system total generation of 5835 MW and loading of 5800 MW were 
considered for the state of load flow analysis. Through the results, it was found that the method is effective. 
To generate the inputs that will feed the artificial intelligence, and about (variables system's bus voltage and 
the angle), multiple-layer intelligence operations were run as the inputs were calculated as the voltages and 
their angles for each bus except for a slack bus. The locations of the phasor measurement units have been 
improved, as in Table 1 about connecting two wind-mills on two buses, the first on Babylon high voltage and 
the second west of Baghdad high voltage. With the presence of air turbines, which were chosen randomly 
location leading to injected of harmonics into the electrical power system, which affects the value of the low 
frequency, and fluctuating it is contained in different values, which leads to instability the Iraqi electrical 
network. To generate a wide range for the system, we review the scenarios by randomly loading the system 
from 100% to 125%. During these operations, a process of outage line from the power network was entered. 
All about those scenarios, the values of our targets such as (MI, DR, and FI) are recorded. 

Several cases for decoration vectors are 1500 cases, 1000 for the training of artificial intelligence 
function, and 500 of them are for testing purposes. After an outage of one line of the power grid lines, the 
data that should be inputted into the artificial intelligence network, which includes both voltages and their 
angles and coming through the phasor measurement units (8 * 4 * 3) 96 are used for each case. Because of 
the number of cases, which is a large number and mathematical process, it will be a matrix with many rows 
and columns, which if we want to analyze will require a long period for this. To reduce the number of 
variables in the required matrix, we use principal component analysis (PCA). To give a 90% reduction with 
an accuracy value of (0.995) and then the data was trained after the reduction using the artificial intelligence 
function (FFNN) to find the target or output of (MI, Fi, Dr) in real-time. The patterns are determined by 
whether it is (local from 0.8 Hz to 2 Hz) or (inter-area from 0.1 Hz to 0.8 Hz) by frequency value. Table 2 
(see Appendix) show the values of each of (MIi, mode rank and mode type) were obtained through in (2) and 
using two methods, each of (PMU-ANN-based method and modal analysis), where they show the difference 
between the results varied in very few proportions. Nine patterns were obtained for each of the three cases. 
The first case 105% (an outage of line KDS-KUT): The loading rate is 105% and the line outage (KDS- 
KUT), as there are 7 cases for one and two who are within the Inter category and for both ways. The second 
case is 115% (an outage of line BAJG-KURK). The loading rate is 115% and the line outage (BAJG- 
KURK), as there are 7 cases for one and two who are in the inter category and for both ways. In a similar 
way for the third case, 125% (an outage of line HADYTH-BAGW) and different results for all cases. From 
Table 3 the values of (Freq.(fi)) and damping ratio (C i), where the values gotten from PMU-NN-based 
method was compared. From the results, it was found that the difference between them was very little and 
with very high accuracy to predict the objectives of them. 

Through the results in Tables 2 and 3 (see Appendix), predictive results were obtained, which 
provides a good and accurate for selecting the type and method of control required in real-time to control the 
damping fluctuation at the frequency of the electric power grid. As is well known to those with specialization 
in electrical power systems, the behavior of the electrical model can change when overloading cases when 
the power system is dynamic. The oscillatory dynamics are obtained in real-time using the proposed method, 
regardless of device size and operating condition. An operator can concurrently figure out the main driver of 
that critical mode or most critical generator for the poorly damped situations, and can, therefore, apply 
effective helpful achievement and use of this approach in real-time also. The MSE (training data) is 0.0156 
and MSE (test data) is 0.0172. 

The MSE is employed as a performance function. From the result, the mean square error is very less 
for both the training and testing phase for the FFNN. The number of samples is 11000 samples. 10000 for 
training and 1000 for testing 


6. CONCLUSION 

In this paper, a new PMU-NN monitoring system for low-frequency oscillation is proposed, the 
linear valuable method is used to the optimal location of PMU to provide full observation of power system 
network buses and reduced the input data to make the convergence of the power system easy. The 
performance of the preparation was confirmed on large-scale power networks. Priority monitoring is given to 
the critical buses; the full information was provided from PMU as input data for FFNN. The artificial neural 
network aims to estimate the low-frequency oscillation, mode type, model index, and damping ratio. From 
the results the proposed method proved can work with high precision in off-line and real-time. The new 
method tested on the Iraqi power system high voltage (400 KV) by using MATLAB package toolboxes and 
compared with another method is modal analysis. 
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Table 2. Different results of mode index and mode rank for model analysis and PMU_NN method for Iraqi 
power system (24 bus) 


Modal analysis 


PMU-ANN-based Method 


Scenario number Loading (%) Mode Mii Mode Type of mode Mii Mode Type of 
No. rank rank mode 
350 (Case I) 105% (Outage of line KDS-KUT) 1 7.263 1 Locally 7.2843 1 Locally 
2 6.7441 3 Locally 6.7654 3 Locally 
3 6.6123 4 Locally 6.6336 4 Locally 
4 6.8278 2 Locally 6.8491 2 Locally 
5 2.9075 7 Locally 2.9288 7 Locally 
6 5.1669 5 Locally 5.1882 5 Locally 
7 3.985 6 Locally 4.0063 6 Locally 
8 1.0301 9 (Inter-Area) 1.0514 9 (Inter-Area) 
9 1.0452 8 (Inter-Area) 1.0665 8 (Inter-Area) 
10 1.0202 10 (Inter-Area) 1.0211 10 (Inter-Area) 
750 (Case II) 115% (Outage of line BAJG- 1 7.3036 1 Locally 7.3236 1 Locally 
KURK) 2 6.5674 3 Locally 6.5874 3 Locally 
| 6.6188 2 Locally 6.6388 2 Locally 
4 4.2785 6 Locally 4.2985 6 Locally 
5 5.163 4 Locally 5.183 4 Locally 
6 3.8456 7 Locally 3.8656 7 Locally 
7 4.8868 5 Locally 4.9068 5 Locally 
8 0.9996 8 (Inter-Area) 1.0196 8 (Inter-Area) 
9 0.9965 9 (Inter-Area) 1.0165 9 (Inter-Area) 
10 0.9887 10 (Inter-Area) 0.9922 10 (Inter-Area) 
1400 (Case IIT) 125% (Outage of line HADYTH- 1 6.9912 1 Locally 7.4012 1 Locally 
BAGW) 2 6.3426 2. Locally 6.7526 2 Locally 
3 6.1593 3 Locally 6.5693 5) Locally 
4 4.0041 7 Locally 44141 7 Locally 
5 5.1446 4 Locally 5.5546 4 Locally 
6 4.8924 5 Locally 5.3024 5 Locally 
7 4.6537 6 Locally 5.0637 6 Locally 
8 1.008 9 (Inter-Area) 1.418 9 (Inter-Area) 
9 1.018 8 (Inter-Area) 1.428 8 (Inter-Area) 
10 1.007 10 (Inter-Area) 0.9936 10 (nter-Area) 


Table 3. Different values for (frequency, damping ratio and DI) of model analysis and PMU_NN method for 
Iraqi power system (24 bus) 


Modal analysis 


PMU-ANN based method 


Scenario Loading Mode Freq. Damping Damping Freq. Damping Damping 
number (%) No. ratio index ratio index 
350 105% 1 1.688 0.071 Poorly 1.6901 0.0742 Poorly 
(Case I) (Outage of line KDS- 2 1.5221 0.0652 Poorly 1.5501 0.0658 Poorly 
KUT) 3 1.3672 0.0723 Poorly 1.3887 0.0732 Poorly 
4 1.3044 0.0765 Poorly 1.3673 0.0761 Poorly 
5 1.2012 0.0881 Poorly 1.2108 0.0877 Poorly 
6 1.0298 0.0863 Poorly 1.0293 0.0861 Poorly 
7 0.9612 0.149 Poorly 0.9379 0.0993 Poorly 
8 0.7594 0.4335 Poorly 0.7332 0.4271 Poorly 
9 0.6591 0.3998 Poorly 0.6461 0.3887 Poorly 
10 0.7903 0.3852 Poorly 0.7863 0.4015 Poorly 
115% 1 1.7725 0.0794 Poorly 1.781 0.0758 Poorly 
750 (Outage of line BAJG- 2 1.572 0.067 Poorly 1.5618 0.0615 Poorly 
(Case II) KURK) 3 1.3943 0.0685 Poorly 1.419 0.0844 Poorly 
4 1.3925 0.0666 Poorly 1.3801 0.0639 Poorly 
5 1.2164 0. 0574 Poorly 1.2915 0.0911 Poorly 
6 1.0023 0.0899 Poorly 1.006 0.0869 Poorly 
ui 0.9316 0.0912 Poorly 0.9458 0.0859 Poorly 
8 0.7903 0.3852 Poorly 0.7863 0.4015 Poorly 
9 0.7133 0.3961 Poorly 0.7001 0.4006 Poorly 
10 0.9692 0.179 Poorly 0.9379 0.0993 Poorly 
1400 125% 1 1.7695 0.0801 Poorly 1.7692 0.0822 Poorly 
(Case IID) (Outage of line 2 1.5615 0.0620 Poorly 1.5626 0.0627 Poorly 
HADYTH-BAGW) 3 1.3937 0.0855 Poorly 1.3872 0.0861 Poorly 
4 1.2833 0.0630 Poorly 1.2825 0.068 Poorly 
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Table 3. Different values for (frequency, damping ratio and DI) of model analysis and PMU_NN method for 


Iraqi power system (24 bus) (continue) 


Modal analysis PMU-ANN based method 
Scenario Loading Mode Freq. Damping Damping Freq. Damping Damping 
number (%) No. ratio index ratio index 

1400 125% 5 1.2111 0.102 Poorly 1.1893 0.1018 Poorly 
(Case III) (Outage of line 6 1.0548 0.0741 Poorly 1.0632 0.0758 Poorly 
HADYTH-BAGW) 7 0.9311 0.0896 Poorly 0.9385 0.0906 Poorly 
8 0.7036 0.0472 Poorly 0.6808 0.0424 Poorly 
9 0.7011 0.0422 Poorly 0.6900 0.0414 Poorly 
10 0.6581 0.3798 Poorly 0.6431 0.3897 Poorly 
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